PROGRAM LIBRARY 


DECUS NO. 

8-549 

TITLE 

POLYNOMIAL LEAST SQUARES FIT 

AUTHOR 

Guy R. Sherwood 

COMPANY 

Edward Hines, Jr. V. A. Hospital 

I Hines, Illinois 

DATE 

I June 28, 1972 

SOURCE LANGUAGE 

PS/8 FORTRAN 


Alt-hough this program has been tested by the contributor, no warranty, express or implied, is made by the contributor. 
Digital Equipment Computer Users Society or Digital Equipment Corporation as to the accuracy or functioning of the 
program or related program material, and no responsibility is assumed by these parties in connection therewith. 










* 


* 








POLYNOMIAL LEAST SQUARES FIT 


DECUS Program Library Write-up 


DECUS NO. 8-549 


Sfu 'STlrtX 13 X dnd tut 

to the main program and both have been stored together*" bd<in linked 

program. The program is stored on rio „i + d together as one core image 

extension. A l? joZll be“j ££tyje! “* “ tha 

V/hen the program has been called and has starts 4 . 

is, "HOW MANY POLYNOMIALS OP DIFFERENT DECREES?"’ J question asked 

o^erent desired polynomials^ JEST ^ger 

different valdes N of I y T for 0 Sch ? x'"" l Thls^Se^™ 0 ? _4L? " ! “ d " H0W MMnr 

fo^at. It is necessary to have botTSe XZXXlXoTlTZ ?. 

there Is IX TvaluJ ST th “ 7*^ ° f thd dat£ ' If 

there can beTly a ma£m£, S 2 XX ^ain. 

each X. The maximum order of the elXll L XXX onST f Y T 
of values vail depend on the Nth derreeTM,,!? T data - size 

experiment to see what is thp w+h ^ * - Thls W1 “ L ~ Va ~ ry and ls bes t to 

rrjssLrjS ----- 

sSr^rvaSfb^ n zxxXs <*r m entered in “ *>»*. 

If no weighting is Ssiref er-IffT ti ’ “”** *" sp “ ed 9* apart. 

First entlr thf XllXttn 7 XlX/ 7^ d 6 is the default 
3 depending). ^ “ elght factor “ d then the Y values (1 through 

After all data has been entereo. the out-v - -htt m . 

the degree of polynomial requested If -'oh n ° W * ^‘ uae ^ii vary on 
entering the data and requestSittat L'“' T™ err °-- happens, try re- 
persists, the Nth dteT^n^LT^e^S^af ^ ^ 


The program will exit to the PS/8 monitor when completed 


vd.tr 


output. 













HOW MANY POLYNOMIALS OF DIFFERENT DEGREES? 03 


HOW MANY DATA POINTS, WHAT DEGREE POLYNOMIAL AND 
HOW MANY DIFFERENT VALUES OF Y FOR EACH X? 

050103 

050201 

040301 

X WEIGHT v 


0 . 0 
0.5 
1.0 
1.5 
2.0 


1.0 2.5 

1.0 4.1 

1.0 6.3 

1*0 7.9 

1*0 11.4 


3*1 2*3 

3*5 4.3 

6.7 6.1 

8.2 g.i 

12.0 11.7 


2 
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NUMBER OF GIVEN DATA POINTS ■ 5 

DEGREE OF POLYNOMIAL = 1 


0 DEGREE COEFFICIENTS 
1 DEGREE COEFFICIENTS 
0 DEGREE COEFFICIENTS 
1 DEGREE COEFFICIENTS 
0 DEGREE COEFFICIENTS 
1 DEGREE COEFFICIENTS 


NUMBER OF GIVEN DATA POINTS 
DEGREE OF POLYNOMIAL = 

0 DEGREE COEFFICIENTS 0, 

1 DEGREE COEFFICIENTS 0, 

2 DEGREE COEFFICIENTS 0, 


NIMBER OF GIVEN DATA POINTS 
DEGREE OF POLYNOMIAL = 

0 DEGREE COEFFICIENT* 0. 

1 DEGREE COEFFICIENTS 0. 

2 DEGREE COEFFICIENTS 0. 

3 DEGREE COEFFICIENT=-0. 


)• 2120000E+01 

• • 4320000E+01 
I.2200000E+01 
•4500000E+01 
•2280000E+01 

• 4320000E+ 01 


= 5 

2 

2577143E+01 
249 1 428E+01 
9 1 428 61E-00 


= 4 

3 

2499999E+01 
1800029E+01 
3599946E+01 
15999 7 6E+01 
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